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Abstract 

The Jackiw-Teitelboim gravity with the matter degrees of freedom 
is considered. The classical model is exactly solvable and its solutions 
describe non-trivial gravitational scattering of matter wave-packets. 
For huge amount of the solutions the scattering space-times are free of 
curvature singularities. However, the quantum corrections to the field 
equations inevitably cause the formation of (thunderbolt) curvature 
singularities, vanishing only in the limit h ^ 0. The singularities cut 
the space-time and disallow propagation to the future. The model is 
inspired by the dimensional reduction of 4-d pure Einstein gravity, 
restricted to the space-times with two commuting space-like Killing 
vectors. The matter degrees of freedom also stem from the 4-d ansatz. 
The measures for the continual integrations are judiciously chosen and 
one loop contributions (including the graviton and the dilaton ones) 
are evaluated. For the number of the matter fields N = 24 we obtain 
even the exact effective action, applying the DDK-procedure. The 
effective action is nonlocal, but the semiclassical equations can be 
solved by using some theory of the Hankel transformations. 







1 Introduction 



The Einstein gravity is undoubtedly beautiful and physically relevant theo- 
retical construction which, at the same time, brought many novel mathemat- 
ical structures into theoretical physics. The notions like the black-hole, the 
horizon, the curvature singularities or the gravitational collapse enriched our 
conceptual world, but they also posed new challenging problems, yet to be 
understood. For instance, what is the fate of matter in the last part of the 
gravitational collapse? Following the overwhelming success of the quantum 
theory in describing the world of subatomic distances, physicists feel that an 
appropriate theory of quantum gravity has to provide the correct solutions 
to the problems and improve our understanding of the phenomena, occuring 
in very strong gravitational fields near the curvature singularities. Unfortu- 
nately, such generally accepted and technically applicable quantum theory of 
gravity does not exist yet. The considerable breakthrough was reached due 
to string theory, where the consistent perturbative S matrix for scattering of 
gravitons and excitations of other fields can be obtained [0. However, the 
dynamics of a gravitational collapse or the status of quantum black holes 
remain nonunderstood. Though there is a common belief that quantum ef- 
fects should smear the singular behaviour of the classical theory, there is no 
sufficiently established quantitative evidence for such a conjecture. 

The string boom had also an indirect, but very important influence on 
the subject of quantum gravity. The models of 1 + 1 dimensional theories of 
gravity coupled to a dilaton field <fi have arisen in string theory ||, ||| . These 
models possess black hole solutions and they motivated Callan, Giddings, 
Harvey and Strominger Q to investigate the U CGHS" model of 2-d dilaton 
gravity coupled to conformal matter. The model is of interest as a "toy" 
model of quantum gravity in two dimensions which contains gravitational 
collapse, black holes, cosmic censorship and Hawking radiation. Moreover, 
the model is very similar to that obtained by the dimensional reduction of the 
spherically symmetric gravitational system in 3 + 1 dimensions, hence, one 
may expect the relevance of the 1 + 1 results to the 3 + 1 physics. Recently, 
many authors have been investigating the quantum dynamics of the black 
holes by using the CGHS model The issues of particular interest are 

the backreaction of the Hawking radiation on the metric and the endpoint 
of the black hole evaporation. The problem is far simpler than the original 
3 + 1 dimensional one and powerful methods of conformal field theory in two 
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dimensions can be used. 

The CGHS model and its variants [15, 16| and also other 2-d dilaton 
gravities have been studied in the literature [H], particularly in the context 
of the noncritical string theory ||18|| . The models can be typically obtained 



by the dimensional reduction of the higher- dimensional pure metric gravi- 
ties. This fact suggests the following CGHS-like scenario for "addressing" 
the four- dimensional quantum problems: one finds the corresponding dilaton 
gravity model in 1+1 dimensions and attempts to quantize it. Though the 
1+1 quantum theory may still be complicated enough to prevent exact solv- 
ability (like CGHS is), usually it is far simpler than its 4-d counterpart. In 
this contribution, we adopt the scenario and address the quantum dynamics 
of the colliding gravitational waves. The fact, that the nonlinear character 
of the Einstein equations results in the formation of curvature singularities 
after collisions of gravitational waves, is known only since seventies [|Tj^, 2U 



and, perhaps, it is less familiar to nonspecialists than the fact that a black 
hole is formed as the consequence of gravitational collapse. However, the 
colliding- waves problem keeps attracting many relativists ||21||-||24|| without 
an interruption, since the discoveries of the first colliding-waves space-times 



by Szekeres [19] and Khan and Penrose [20]. The problem of main interest for 
us will constitute in the following: What is the quantum status of the those 
scattering space-times? As we have mentioned above, one usually expects 
that the curvature singularities should be smeared by the effects of quantum 
fluctuations. Our quantitative analysis will show, however, the surprising 
result: the quantum curvature singularities are even worse than the classical 
ones and even classicaly nonsingular spacetimes are destabilized by quantum 
curvature singularities. 

Apart from the physical questions which our analysis will try to answer, 
the model to be considered in this paper is of interest also for some more 



theoretical reasons. Indeed, while in CGHS and related theories [F7, 1 
the matter degrees of freedom are added by hand after the dimensional re- 
duction, in our model the matter degrees of freedom also stem from the 
four- dimensional theory. This fact should even increase the relevance of our 
results for the 4-d case. There is another pleasant thing, namely, not only 
the matter loops, but also one loop dilaton and graviton contributions can 
be evaluated, yielding the one loop effective action. We need not to perform 
1/N expansion, but we can take h as the natural parameter of expansion. 
Moreover, for the critical number of matter fields (N=24) our result will be 
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nonperturbative and exact. But the good news is not exhausted by that, it 
turns out, moreover, that the semi classical equations can be solved and the 
behaviour of curvature singularities is under control. 

The plan of the paper is as follows. In section 2 we introduce the 2-d 
matter-dilaton model motivated by the dimensional reduction of the 3 + 1 
dimensional system with two commuting space-like Killing vectors. Then we 
find the classical equations of motion in the conformal gauge. The dilaton 
field turns out to satisfy the standard d'Alambert wave equation, hence, we 
introduce a sort of the "light cone" gauge. In this gauge the matter fields obey 



the Gowdy cylindrical wave equation ||25|| , the general solution of which can 



be given by the decomposition into the Fourier-Bessel and Fourier-Neumann 
modes. The corresponding metric we find explicitly by integrating the re- 
maining equations. We show that the Neumann modes cause the formation 
of the (classical) curvature singularities which close the space-time to the fu- 
ture, while the appropriate superpositions of the Bessel modes describe the 
collisions of the wave packets travelling against each other with the velocity 
of light. The corresponding spacetimes are everywhere regular with the out 
region in which the scattered wave packets travel to the opposite space infini- 
ties. In section 3 we discuss the quantization of the model. We choose the 
standard Polyakov measure for the functional integration over the metrics 
and the reparametrization invariant measures for the dilaton and the mat- 
ter fields integration. Then we compute the one loop effective action . The 
effective action is nonlocal even in the conformal gauge, due to the presence 
of the direct matter-dilaton coupling in the action. The one loop effective 
field equations are localized by going to the dilaton "light-cone" gauge. The 
renormalization requires a purely dilatonic counterterm, the contribution of 
which makes finite one infinite constant in the semiclassical field equations. 
The actual computation requires the knowledge of the functional derivatives 
of the determinant of the (Gowdy) wave operator with respect to the dilaton 
and metric. They are evaluated by using the heat kernel regularization and 
some theory of the Hankel transformations in the Appendix. In section 4 we 
solve the semiclassical field equations. We perform the detailed analysis of 
the scalar curvature of the space-times, obtained by solving the semiclassi- 
cal equations. We show that the contribution of the quantum fluctuations 
to the effective action inevitably generates the curvature singularities in the 
semiclassical spacetimes. These singularities may disappear only in the limit 
h — > 0, thus indicating, that the classical regular scattering space-times are 
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in fact unstable from the quantum point of view. We end up with short 
conclusions and outlook. 



2 The model and its classical dynamics 
2.1 The dimensional reduction 

The form of the four-dimensional metric describing the collisions of collinear 
gravitational waves is given by |2^] 



ds 2 = -2<f)-h^dudv + <f){e V ^ Q dx 2 + e-^ Q dy 2 ), (1) 

where the metric functions /i, <fi and Q are invariant on the (x, y) plane of 
symmetry. The 4-d vacuum Einstein equations for the metric ([!]) consist of 
the constraints 



- 4>uu + Hu4>u = K(J)Ql (2) 

-cj) vv + fi v (f) v = K(f>Ql (3) 
and the evolution equations 

0™ = (4) 

{<t>Qu)v + (<PQv) u = o (5) 

It is not difficult to demonstrate, that the same set of the constraints and 
the evolution equations follows from the following 2-d action 

S = Yk I ^^9^ R ~ KQ^daQdpQ), (7) 
after fixing the conformal gauge 

ds 2 = -2e M dMdw. (8) 
The action (|7]) can be interpreted as the Jackiw-Teitelboim gravity [26] where 



the cosmological constant is replaced with the kinetic term of the matter. 
The matter is coupled to the dilaton and possesses all dynamical degrees of 
freedom of the theory. 
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2.2 The solutions of the field equations 



In what follows, we shall consider the model (0) with an arbitrary number of 
matter fields. The classical dynamics does not change dramatically, but the 
properties of the quantum theory will depend on that number. The action 
reads 

S = h I d ^^^ R - K 9 aP d a QidpQ j ) (9) 

and the constraints and the evolution equations in the conformal gauge get 
obviously modified 



vv 



+ fx u( f> u = K(f>QiQi (10) 
+ /i v (j) v = K<f>QlQi (11) 

<puv = (12) 

t)v + Wi)u = (13) 

- n uv = KQiQl (14) 
The general solution of (|T2"D reads 

<j> = f(u)+g(v) (15) 

If <f> v = (or <j> u = 0), then from QT1]) (or (0)) it follows that Ql = 
(Q 3 U = 0) anc l from (0) \i uv = 0. Since the scalar curvature R is given by 

R = 2e~y uv (16) 

and in two dimension the curvature tensor reads 

Rafl-fS = 2(9a^9l38 — 9aSgi3y)R (17) 

we may conclude that arbitrary functions 4>{u) and Q->{u) (or cf>(v) and Q-*(v)) 
are solutions of the field equations and the corresponding space-time is flat. 
Such solutions obviously describe the matter excitations propagating in one 
direction with the velocity of light. 

If neither <p v = nor <j) u = we can (at least locally) perform the confor- 
mal transformation 
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U = y/2f(u),V = y/2g(v) (18) 

Such a transformation is obviously the symmetry transformation of the set 
of the field equations (|T0|-|T^) , hence we may fix this residual symmetry by 
the claim 

U + V 

#=-^— < (19) 

We may call this gauge fixing the "dilaton" gauge. In the gauge the field 
equations become 



t*u = <U + V)Q\jQ\j (20) 

»v = k{U + V)Q 3 V Q ] V (21) 

Qit + - t Q\ - QL = o (22) 

- fi uv = kQIjQ j v (23) 



where 



We observe that the linear equation for the matter fields Qi does not contain 
the metric function \i. We may call this equation by the name of Gowdy, 
who studied cosmological models with plane symmetry |25| , governed locally 
by (|2"2"D. The general solution of the Gowdy equation (which tend to zero at 



the spatial infinities) is given by the following mode expansion [24, 2S . 



Q j (t,o-)= dujRe[A J (uj)J (ujt)e- iuJ ' 7 + B^(uj)N (ujt)e- iuJ ' 7 ], (25) 
J o 

where A 7 and are (complex) distributions ensuring the proper behaviour 
at the space infinity and Jo and iVo are the Bessel and Neumann functions 
of zero order, respectively. This mode expansion can be easily found by 
using the Fourier tranformation in the variable a in (^2|). The resulting 



6 



ordinary differential equation in the variable t is then the Bessel equation. 
We should note, at this place, that in higher dimensions some additional (so 
called "solitonic") terms are considered at the r.h.s. of (^Sj). They do not 
vanish at the space infinity and in the limit of the weak gravitational coupling 
k — ► those solutions diverge and do not approach the non-interacting matter 
solutions |30| . We shall not consider this "solitonic" sector in this paper and 
prescribe the boundary conditions, mentioned above. 

It remains to solve the equations (pO|),(f2lD and (|23|), which determine the 
metric function //. Combining (|20D with ([H]), we obtain 



fi a = 2ntQ{Qi 



(26) 



fit = Kt(QiQi + QiQi) (27) 
We use the fact [27] that for F and G, fulfilling the Bessel equations 



it holds 



A 2 F 



dF 



x 



dx 2 ~^ dx 



+ (AV - n 2 )F 







dx z dx 



dG 



. dF . 



(28) 



dx(X 2 - u 2 )xF G = x(F^-G- 

1 dx dx 

This formula enables us to integrate the products of the Bessel functions. 
The result of the integration gives the explicit form of the metric functionP] 



k I d(Jiduj2UJiUJ2 t Re 
h 



1 



G{{u 1 t)G J {u 2 t)e- i ^+^' 7 



1 1 



2 u)\ — UO2 



(G{* (w 1 t)G J (cj 2 t)e i{wi -" 2)CT - G{* {^GHu^e-' 1 ^-^)] , (29) 



1 It appears, that this result is new even from the point of view of 4-d theory of colliding 
waves. 
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where 

G J 0{1) (ujt) = A j (u;)J 0{1) (ut) + Bi{u)N Q{x) {ut) (30) 

We note, that the classical equations ( |i~2| - [r4l ) turn out to be "iteratively" 
linear. Indeed, solving the linear Eq . fll2|) and inserting its solution into 
Eq. (|i~3|) , we get again the linear equation. After solving it, we insert Qi 
into Eq. flUP and get the linear equation for fi. Such a structure of the 
equations gives the classical integrability and will be also important later for 
the quantization. 



2.3 Curvature singularities and the global structure 

We start our analysis of the curvature singularities with the formula for the 
scalar curvature. Following from Eqs. (|T6| ) and (^3|), we have 

R = Ke -»(-QiQi + QiQi). (31) 

Near t — > + we have 

J (t)~l-^ + ... (32) 
N Q (t)^(l- t l)\nt + --- (33) 

hence, 



1 r°° E j 
Ql ~ -( / duuRelB 1 (u)e~ 1LVr7 }) + bounded = — + bounded, (34) 
t Jo t 

rod 

Q{ ~ lnt ( / duuRe[-iB j (u)e-' luJa ]) +bounded = H j lnt + bounded, (35) 
Jo 

(i ~ k In t E j E j + bounded. (36) 
Inserting Eqs.(|34"D,(|3"5D and (^) into Eq. (|3l|) , we have 

R ~ r KEJEJ [ — + WH j (\nt) 2 + ■■■). (37) 

We conclude, that the regularity of the (classical) space-times requires both 
Ei and if J to be equal zero, or, equivalently, 
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B j (uo) = 0. 

Consider now (regular) space-times, given by 



(38) 



A 4 (uj,ujo,p) = lajle^ , 



-e 2 p 
Ap 



BHuj) = O, 



(39) 



where <pj is real and p and uq are real positive parameters. Note that for 
p^O 



A^[uj, u , p) — > \a,j\ei(j)j— ^i=5{uj — Uq). 

v2 



From Eq.(P5|) for = we obtain 



Q J 



dcjRe 



\aj\ern 



-e 2 p Jo(ujt)e 
Ap 



(40) 



(41) 



From the well-known formula for the asymptotic behaviour of the Bessel 
functions for t — > ±oo p7 



we have for t — > ±oo 



2 , 7T . 

cos {ut =F -r) + 



HUJ \t 



(42) 



li/e-^*-' 7 ) 2 cos [uo{t -a) + 0^ =p 



7T n 



+ e - §(*+-) cos[a;o(t 



(43) 



Now the physical interpretation of this solution is obvious. At t — > — oo 
we have two wave packets propagating against each other by the velocity of 
light; at t — > oo the two scattered packets propagate apart from each other 
with the gained phase shift, indicated in Eq. (f4"3"D . Because Joiut) and its 
derivatives are bounded functions [27|, we may use the Riemann-Lebesgue 
lemma and from Eq.(41) conclude that for t = const and a — > ±oo, Qi 
and all its derivatives with respect to t and a vanish. Hence, by using the 
constraints (|2~6D and and the formula (|3~1~D for the scalar curvature, we 
observe that the space-time is flat in this limit. For the cases a = const, 
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Figure 1: The scalar curvature of the regular classical space-times. 



t — > ±00; t + cr=const, t — cr — ► ±00 and t — cr=const, t + a — > ±00, we use the 
asymptotic expression fl43|) , the constraints (^), (|27|) and the formula (j3ll) 
to arrive at the same conclusion. Therefore, for the "wave-packet" choice 
([39D the corresponding space-time is asymptotically flat (see Figure 1), it 
has the same topology as the two-dimensional Minkowski space-time and is 
free from curvature singularities. We shall not need the explicit form of the 
metric, which, nevertheless, can be obtained by performing the integral ( p9|) 
with the choice (|39|). We should end up the classical analysis with some 
important comments. 

First of all, it does not seem unexpected that for the collisions of the 
localized wave packets travelling with the velocity of light, the space-time 
is asymptotically flat in the space-like and time-like directions. What looks 
more surprising is the fact that the same is true for the null infinities. The 
reason is simple: in two dimensions a single propagating wave does not curve 
the space-time (cf. the analysis after (Eq.|P7D). In higher dimensions this 
is not true |29|, [30|1 , but in that case the curvature is given by the shape of 
the wave front in the transverse directions. Since there are no transverse 
directions in two dimensions, our result could be anticipated. 
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The second comment is closely related to the first one. It concerns the 
regularity of the initial data. In the higher dimensional case the following 
problem was studied [^3], |SD[. If initially regular gravitational waves interact, 
will a curvature singularity be formed? The criterion for the regularity of the 
incoming data can be naturally formulated: one requires the boundedness of 
the amplitude of the wave,which itself is defined by means of the components 
of the Riemann tensor of the corresponding metric in the so called parallelly 
propagated orthonormal frame p3[ |3(J. However, in the two-dimensional 
case, incoming waves do not curve the spacetime and this criterion of reg- 
ularity fails. But certainly we should not consider all incoming waves as 
regular (with bounded amplitude). One might require that the scalar fields 
Qi itself should be bounded, but, on the other hand, we could rescale it by 
an arbitrary function of another scalar field, the dilaton 0, and we would 
get the classically equivalent dynamical theory with the different condition 
of incoming regularity. Fortunately, it appears to be a natural candidate for 
the amplitude of the wave. The matter part of the action @ suggests the 
following inner product on the space of fields Q (see also Eq.(|5Tj|)) 

(Qi,Q 2 ) = / d 2 £y=#Q!Q 2 . (44) 
Hence, our condition of the incoming regularity reads 

lim (j)Q j Q j = finite (45) 
In the case of the wave-packets (j39|) we get 

^lim <j)Q j Q j = |ai| 2 e- 2py2 cos 2 [V2u; 2 V - fa T \]= finite (46) 
and similarly for V — > ±oo. 



3 The quantization 

3.1 The functional measures and the effective action 

Define the generating functional of the model W[J g , J^, Jj] by 

W[J g , J„ Jj] = J Vol(w/) (47) 
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exp {i J d 2 £[^^(R - Kg a/3 d a Q j d p Q j ) + ^J 9 R + + JjQ 3 ]} 

where J g is a scalar source, and Jj are scalar densities and Vol(Diff) 
is the volume of the group of diffeomorphisms. We define the functional 
measures by the following norms 

I N I 2 = / d 2 £V^ a V*<W<^> (48) 

INH 2 = / d 2 ^H 2 , (49) 

II^H 2 = / d^v^^W- (50) 

Eq.(|4Tf) defines the usual de Witt-Polyakov norm |JT] and Eq.(^) gives the 
standard reparametrization invariant measure for a scalar field. The norm 
( |5T)D is given by the form of the matter part of the classical action, much in 
the same way as the norm for the quantization of the standard non-linear 
er-model with the coordinates X A of the target and a metric Hab(X), i.e. 

[0 

\\6X^\\ 2 = J d 2 £V=gH AB (X(0)5X A (05X B (0. (51) 
We return to Eq. (|47|) and we fix the conformal gauge 

ds 2 = -2^dudv. (52) 

By using the standard Faddeev-Popov procedure, we obtain 

r 26 r 1 

W[J] = j D^D^D^Q* exp {i— J d 2 £-/i«9V} (53) 

x exp {- f d 2 £[^(-«9 2 /i - KdQ j dQi) - J g d 2 fi + J^> + JjQj}}, (54) 

where d 2 is the Minkowski d'Alambertian and dQWQi means the Minkowski 
metric scalar product. The Weyl anomaly term comes from the Faddeev- 
Popov determinant and the measure D^n is given by the following norm 

|M a = |dVV- (55) 
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We standardly suppose that the exponential term in the Weyl anomaly is 
eventually (after including all contributions) cancelled by tuning of the 2-d 
cosmological constant counter-term. 

Define now the effective action T by the following prescription 

y In W[J) = Z[J) ee I> c , <j> e , Q{] - J g d 2 fi c + J^ c + JjQl (56) 

where 

,9V ee- — , C =— , Q J c = JT- (57) 

OJg OJtf) OJj 

We wish to compute the one loop effective action IV In order to do that, we 
have first to determine Z\ (the generating functional for the connected Green 
functions) from Q53D and then to perform the Legendre transformation fl5E|) 
and (|57|). We note, that the dependences of the measures on the fields \i and 
cf) are of the order 0(h) with respect to the classical action in the exponent. 
Hence, the loop diagrams with the vertices coming from the measures will 
be of the order 0(h 2 ) and may be neglected in the one loop approximation. 
Therefore, we may write 



f 26 f 

W S emici.[J) = J D^j/iD^cpD^j^jQ 3 exp{i— J d 2 £ fijd 2 fij} 

exp J d 2 e[^(-«9 2 /i - KdQ j dQ j ) - J g d 2 fx + + J 3 Q j ]}, (58) 

where //j and <pj are the saddle point values of the exponent, given by the 
equations 

- <9 2 /ij - KdQ j dQ j + 2kJ^ = 0, (59) 

0j + 2nJ g = 0, (60) 

d((pjdQ>j) + J, = 0. (61) 

We observe that in the one loop approximation, we can consider the measures 
to be independent on the field integration variables (but, of course, depen- 
dent on the Schwinger currents). Now we evaluate the integral (|58f) . In the 
(second) exponent, there stands the quadratic form in the variables <ft an d 
jj,. Moreover, the norms defining the measures has the same /xj-dependence., 
i.e. 
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|M| 2 = J d 2 £e^ 2 ; \m 2 = J dV J ^ 2 - (62) 

Therefore, we can easily perform the Gaussian integration over fi and <fi with 
the result 

r 24 r 

W semicL [J} = J D^j^jQ 3 exp{i— J d 2 £ /ij<9 2 /ij} 

exp | d 2 ^{-2Kj g J^ + nJ g dQ 3 dQ 3 + J,Q 3 )}. (63) 

The integration over Q 3 is again Gaussian, hence we obtain the closed ex- 
pression for the semiclassical generating functional 

W semicL [J g , J+, J,} = det~ N/2 [ - d e -^-^—d(J g d)} 

lb ZsHjU g 

where we used the definition ( [5"0|) of the measure D^Q 3 and Eq. (|60D . The 
Legendre transformation (^) and (|57D can be easily performed and we obtain 
the following expression for the one loop effective action 

r^^Qi) = f d 2 A-d 2 /i c - KdQidQi) + n^- [ d 2 ^ c d 2 f, c 

J Ik 9o7r J 

+i/Andet[ - T^-e'^^-d^cd)} + 0(h 2 ). (65) 

We recognize the classical action and two quantum corrections. The first one 
is the Weyl anomaly, while the second one is the nonlocal term depending 
on /i c and C . The expression (|65| ) can also be written in the manifestly 
covariant way, i.e. 

r x (^, C , Q 3 C ) = J d 2 £^ ^-{R c - KgfdaQidpQi) 

24 r 1 

96tt J y^g c 

+iftylndet[ - ^-—L-dai^^gfdfs)} + 0{h 2 ). (66) 
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3.2 The case N = 24 



For N = 24, the quadratic term of the Weyl anomaly vanishes. The dilaton 
gravities usually get simplified and more precise results can be obtained in 



that case [Q. This happens also in our model. We show that the semiclassi- 
cal effective action Ti fl65|) is the exact quantum effective action of the theory 
for N = 24. We use the DDK-approach P4| |5o]| to establish this result. The 
dependences of the functional measures on the field /x read 

ZV = (Do») exp / d 2 ^fid 2 fi], (67) 

Dfi(f> = {D cf>) exp J d 2 {~f,d 2 t,), (68) 

D^Q 3 = (PwQ?) (69) 
X exp [ - -±- J d 2 £{^V + ^[2d 2 In |0| + (d In |0|) 2 ]}]. 



The formulas (|6"?D and (|68]) are fairly standard |36|, |57| , however, the relation 
( |69|) deserves some comment. Indeed, it can be explicitly derived by com- 
puting the Jacobian, which relates both measures, with some regularization 
procedure. We shall use the heat kernel regularization and use the defining 
formula ( p0|) to write 

D M Q = Do^QVfetL, (70) 
where L is the diagonal operator 

£(6,6)=e^ l} 5(£i,6). (71) 

Note that the 5-function 5(^1,^2) is to be understood in the sence of the 
scalar product floTf ) with /1 = 0, i.e. 

= -7^*& -6). (72) 

Clearly 

51ndet L = 5Trln L = J d 2 £0(O<K£, OMO, (73) 
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where is the meaningless quantity.As in |36], |37|, we replace it by the 

heat kernel of the covariant Laplacian, but in our case with respect to the 
scalar product (150), i.e. 



8 e {£, = -jr < £| exp { - e[ - l-^d a {yf=gif^dp)\}\^ >- 



']}• (74) 



We obtained the last equality in the conformal gauge, by combining the 
formulas ( |115| ) and ( |121| ) of the Appendix. Now we insert ( |74|) into ( |73|) and 
in a straightforward way we arrive at the formula (169|) . 

We may also check the validity of the formula fl6"P| ) for a particular inte- 
grand. Indeed, let us compute the integral 



D^Qexp | d^Qd^Q] = def- 1/2 [ 



27i 



-9(09)]. (75) 



We have (see Appendix Eqs. (|122| , |123D ) 
def 1 ^ 1 

x exp { - -±- J d^i^/i + ^[29 2 In |0| + (9 In 



2ft 



^(09)]=det- 1 / 2 [--^(^)] 



1}}- 



i 1 

2h'< 



Because 



(76) 



(77) 



Eqs.([75D, ( |76D and (|77|) obviously match with the formula (p9|) . 

After this digression, we now compute the effective action for the case 
N = 24. We use the defining formula (|54]) for the generating functional 
in the conformal gauge and insert the field dependences of the measures 



det-n-^W)] 



A^Qexp[— / d 2 £Q«9(0«9)Q], 



(|3), (HD and (H) in it. We obtain 



W[J] = D fi D Q <j> D ^Q j exp{ 



24i 
32^ 



d 2 £/i[2«9 2 ln|0 c | + («91n|0 c |) 2 ]} 



exp | d 2 £[^;(-<9V - KdQ j dQi) - J g d 2 fi + + J.Qj}} 
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(78) 

The integration over Qi is Gaussian and over /i gives the 5-function, therefore 

W[J) = [ D <p5(d 2 <P + 2Kd 2 J g + ^[2d 2 ln\<f ) \ + (d\n\<f ) \) 2 ]) 
J 2tx 

xdet -[ _ lA am] exp _ j , i m = 



det " 12 i-2k^) 9 <™ 



x«p[i/d'«w,)- ^^a)^)]. P9) 

where the dependence of 0(J S ) on J 9 is dictated by the 5-function in (|79|). We 
stress that the formula ([79|) gives the exact generating functional. Performing 
the Legendre transformation and (|57|) we obtain the exact effective action 



r(/i c , C , q>) = | d 2 ^(-«9Vc - Kflgjag*) 

-i- / d 2 £/i c [2«9 2 ln|0 c | + (<91n|0 c |) 2 ] + m% lndet[- ~d{<j, e d)\. (80) 

47T J ZAi (B r 



Comparing the result with Eq.(^), we conclude that for iV = 24 the semi- 
classical approximation is, in fact, exact. 



4 Quantum curvature singularities 

4.1 The semiclassical field equations 

We obtain the semiclassical field equations by varying the one loop effective 
action 1^ with respect to the classical fields fi c , C and Q 3 C . We have 

- + h w^ + h Yw c lndet[ - ^ c j c d{m] = °' (81) 

d(<f> c dQi) = 0, (82) 
- d 2 [i c - KdQidQi + HiKN— lndetf - -U" Mc — dU c d)\ = 0. (83) 

OQr 2fl r 
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The solutions of the semiclassical equations have expansions 



He = Hc,o + hfic,i + O (h 2 ), (84) 

c = c ,o + ft0c,i + O(ft 2 ), (85) 
Qi = Qi >0 + hQi tl + O(n 2 ). (86) 

Because we know just the first loop corrections to the effective action, the 
0(h 2 )-terms in the field expansions (|8jj) ,(85) and ([86]) are irrelevant in the 



one loop approximation. Our next task will consist of the determination of 
Hc,i, 4>c,\ an d Qii from the semiclassical equations (|81"|),(|82|) and (|83|), when 
Hc,o, c ,o and Q 3 c is a given classical solution. Since the "lndet" terms in the 
one loop field equations are already of the order 0(h), it is enough to compute 
the functional derivatives of lndet at the classical solution /i c and c ,o(= t). 
The actual calculation requires the knowledge of the heat kernels of elliptic 
operators, some theory of the Hankel transformations and some integration 
of the Bessel functions. The details are presented in the Appendix here we 
list only the final result 

^" ldet l4 e ""i a( ^l = ^ +I ^[2#ln W +(31nW) 2 ]. (87) 



i_ e -"-i-9(0 o a)]U =1 = (88) 

i i fit 2 
— \d 2 a c - d(iicdln \t\)] .In (— ), 

where Q is a finite constant and e — > is the ultraviolet cut-off. The following 
counter-term is needed for the cancelation of the UV divergence 

C cmnt . ~ h^ g af3 <p- 2 d a <pdp<p. (89) 

This is the vertex of the new type, as it is usual in the renormalization of 
2-d field theories. Its appropriate tuning replaces the product efl in the 
semiclassical equations by some finite constant, but it cannot remove the 
logarithmical dependence on t 2 . 



———lndet [ — 
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Inserting the evaluated functional derivatives (^7]) and fl88| ) and the h- 
expansion (jS4|) , (|85]) and (|86D of the fields into the semiclassical equations 
flgT| ) , (^) and ( PB| ) we obtain the equations for /i Cj i,0 Ci i and Q 3 cl 

1 9 24 -W 9 N 

- + —r^—d 2 ^ - -— = 0, 90 

d(<t>c,idC&,o) + d(tdQi il ) = 0, (91) 

(-5 2 /i c ,i - 2k8Q; 1 x dQ\ 1 Q ) = ^(7<9 2 /i c ,o + i^t^co - i^c,o - i In— ^— ). 

87r i t z const 

(92) 

This system of equations can be solved in a similar way as the classical 
system (|T2|) , (|13|) and ( |i4|) was solved. Indeed, because we know the Green 
function of the Minkowski d'Alambertian, we find from Eq.(^) the general 
form of Cj i, by adding arbitrary solution of the homogeneous equation to one 
particular solution of the full equation. Inserting c l into (pTj) , we obtain 
the linear (inhomogeneous) Gowdy equation for Q J cl . Since we know the 
eigenvalues and the eigenfunctions of the Gowdy operator, we know also its 
Green function and, eventaully, the general form of Qii- Finally, putting 
Q J c i into (|9"2"D, we obtain the linear inhomogeneous d'Alambert equation for 
/i Cj i, the general solution of which can be easily found. We conclude that 
our semiclassical equations (|90|) , (pl|) and (^) are exactly solvable. For 
our purposes, there is no need to write down the explicite (and somewhat 
cumbersome) formulas. Instead of that we shall concentrate on the behaviour 
of the general solution near t ~ 0. We shall show, somewhat surprisingly, 
that even when we consider a regular classical solution /i Cj0 , c ,o and Q 3 C>0 , the 
corresponding solution /i Ci i,0 Cj i and Q 3 cl possesses necessarily the curvature 
singularity at t = 0. Such space-times are therefore classically regular but the 
quantum fluctuations induce the scalar curvature singularity, proportional to 
h. Hence, the quantum effects not only do not smear the classical curvature 
singularities, they even destabilize the regular space-times! We present the 
corresponding quantitative analysis in the next subsection. 



4.2 The scalar curvature of the semiclassical space- 
times 

Let us study the behaviour of the scalar curvature R near t ~ for the 
space-times which solve the semiclassical field equations. In this subsection 
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we omit the index "c" of the fields /i c , C and Q{. We choose such classical 
metric field /j,q, that the classical space-time is nonsingular. From Eqs. (|25|) 
and (^) for B j = 0, it follows for t ~ 

/i ~t 2 /(^) + ---, (93) 
Qo = 2^)+t 2 /^) + ---, (94) 

where f(a),g(a) and fo(cr) are functions, the concrete form of which is not 
relevant for our purposes and dots denote the subleading terms, also irrele- 
vant for our analysis of the curvature singularities. Hence, from the Eq. (|90|) 
we find the behaviour of 0i near t ~ 

Nk 

' In ]t] + const f(a)t 2 + F(U) + G(V) + ---. (95) 
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The functions F(U) and G(V) cannot be specified from this equation, how- 
ever, we may change our "dilaton" gauge condition ([19]) by the prescription 

t + h[F(U) + G(V)]^>L (96) 

The one loop effective action (^) is invariant under this transformation, 
hence, the semiclassical equations (|8T|),(|82|) and ( |83| ) remain unchanged. 
Moreover, the classical solution at which the functional derivatives of the 
determinants are to be evaluated change just by the terms of order h. This 
effect, of course, remain unseen in the one loop approximation, because the 
"lndet" terms are already of the first order in ^-expansion. Thus, all sub- 
sequent analysis goes through and we can omit the F{U) + G{V) term and 
write without a loss of generality 



In |t| + const t 2 + ■■■ (97) 



167T 

Now we insert <p% into Eq. (|9~lD . We obtain 

, d 2 Id d 2 H(a), ,, Nnhia) , , , 

r ~dC- ~ tdi + = -T 1(1 " + b0Unded ' (98) 

where H(a) is some function of irrelevant shape. It is easy to determine the 
behaviour of a particular solution of (p8|) near t ~ 0. It is given by 



Qi, par = ~H{a)t In |t| + [2H(a) + N *f(°> ]t + ... (99) 

07T 
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A general solution of the Q\ equation near t ~ is then given by Eq.(|99|) 
plus an arbitrary solution ( p5|) of the homogeneous Gowdy equation. From 
and (0) then easily follows that 

dQidQl ~ C{a) In |t| + bounded. (100) 

The function C(a) vanishes, if the Neumann modes are absent in the "ho- 
mogeneous" part of Qi. Inserting (|100|) into the remaining semiclassical 
evolution equation (P2]), we obtain 



- <9 2 /ii = 2kC(<t) In \t\ - + I l n + bounded, (101) 

8n t t A const 

hence 

Nk Tit 2 

ui = In + p(U) + u(V) + bounded. (102) 

167rt const 

In the classical case, the arbitrary integration functions p(U) and v(V) 
are determined from the constraints (20) and ([H]). In the semiclassical case 
they have to be determined from the given boundary conditions. The situa- 
tion is fully analogous to that occuring in the CGHS model f|], where the 
semiclassical contribution to the constraints come from the Polyakov nonlocal 
action 

s ? = +h ^— / ^V^gR(^d a V^gg a %r 1 R, (103) 

9o7r J \ — a 



9Qtt J V~9 
We can make the generaly covariant action (|103|) local at the cost of intro- 
ducing a new auxiliary field Z (in a similar but not identical way as in |14| ). 
It reads 

94 _ Af r 1 

S P = h^^- J & 2 tV=g[^g aP d a Zd p z + RZ). (104) 

The contribution to the constraints are then obtained from Sg uu and 5g vv 
variations of the action, in the conformal gauge. We have 



1 5S P % 2i-Nl 



r=g 5g uu 48vr v 2 

1 5S P ^ 24 - N , 1 



h—- (-Z V Z V - Z vv + p v Z v ), (106) 



y/=g5g vv 48tt v 2 
We get rid of the auxiliary field Z using the equation of motion 

d 2 (p + Z) = 0, (107) 
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hence 
and 



fi = -Z + p + {u) +fi~{v) (108) 

5S P ,24 — N 1 . I . . . , 

n—Z [Pun ~ 7>HuHu) + T^(u), (109) 

457T 2 

24 - iV 1 

fi—rz {Hvv - ■zfJ-vfJ'v) + T (v), (110) 



^g5g uu 48vr ^ 2' 

1 8Sp _ 1 
v/=p5<T ~ ~48^r (Aiw ~ 2' 

The functions T + {u) and T~(t>) are undetermined, because p + {u) and p~{v) 
are not known. 

In our present model the Polyakov nonlocal action is the part of the one 
loop semiclassical action (c.f.fl66D). Therefore, the unknown functions p(U) 
and v(V) can be specified only by fixing the boundary conditions. The 5g uu 
and 8g vv variations of the remaining lndet part of the effective action of our 
model cannot influence this conclusion and we shall not consider them. 

Finally we are ready to write down the scalar curvature of the semiclas- 



sical space-times. It reads (see Eq. fllq) ) 

R(k) = - e -^d 2 n(h) = R(0) + h(-R(0)iMi - e-^d 2 ^) + 0{fi 2 ) = 

Nk fit 2 

= «0» + M«(o)— to— -jko) W io + ,oo)} 

Nk 1 fit 2 

+h{ - e-"2KC(*) In \t\ + e-^ — (f(a) + ^^)} 

+bounded + 0(ft 2 ). (Ill) 

Clearly, whatever the functions p(U) and u(V) may be, the semiclassical 
space-time is obviously singular. The singularity occurs at t — and all 
timelike observers will run into it - which is referred to as a "thunderbolt" in 



Ref . ||TT[j . We arrived at the remarkable conclusion, that while at the classical 
level there existed the nonsingular spacetimes, at the semiclassical level all 
space-times are necessarily singular. From Eq.fl37|) and (|111|) we also learn 



that the singular behaviour of classical and quantum curvature is different, 
hence no cancelation of a classical curvature singularity due to quantum 
effects may occurfj If the classical spacetime is regular, then the formula 

2 I am grateful to R.Jackiw for a comment on this point. 
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( |1 1 1\ ) says, that the corresponding semiclassical space-time is plagued by the 
curvature singularity proportional to %. Schematically 



R = regular + % singular ■ ■ ■ . (H2) 

We conclude, that the quantum effects destabilize the classical space-times 
and lead to even more severe curvature singularities then the classical dy- 
namics does. There remain only one possibility to avoid this conclusion in 
the framework of the present model, which may seem quite unnatural, how- 
ever. It consists in introducing by hand into the effective action several finite 
counter-terms of new type, which would be fine-tuned so that to cancel the 
divergent terms in ( |101|) . But also keeping this possibility in mind we may 
conclude that the quantum instabilities in our model are generic. 

5 Conclusions and outlook 

We attempted to give a detailed description of the classical and quantum 
dynamics of the Jackiw-Teitelboim gravity with the cosmological constant 
replaced by the kinetic term of matter fields. We showed that the clas- 
sical solutions of the model have natural physical interpretation, namely 
they describe the collisions of the wave packets of matter. For huge class 
of such solutions the corresponding classical space-times are topologically 
trivial, asymptotically flat and free of curvature singularities. Then we com- 
puted the semiclassical effective action of the model; for the case N = 24 
we, in fact, got the exact expression. The effective field equations turned out 
to be manageable from the technical point of view. We have solved them 
and provided a simple analysis of the semiclassical solutions near t ~ O.The 
surprising result followed: the scalar curvature acquires the quantum cor- 
rection which is necessarily singular. Hence, the quantum fluctuations do 
not smear classical curvature singularities, in fact they do just the opposite 
thing, they plague the regular classical space-times with quantum curvature 
singularities. Because for N = 24 we obtained the result starting from the 
exact effective action, our conclusion does not seem to be an artefact of the 
semiclassical approximation. 

We believe that the model which we investigated is also interesting from 
the field theoretical point of view. At the classical level it is completely 
integrable and iteratively linear in the sense of subsection (2.2). This kind of 
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"linearity" played the decisive role in the evaluation of the continual integral, 
in a similar way as it was reported recently in the context of 2+1 Chern- 
Simmons theory The fact that for iV = 24 that computation gives 



the exact result suggests an existence of a deeper algebraic structure in the 
modefj Moreover, it turns out that the model possesses an unexpected and 
interesting geometric structure. Indeed, the action (|7|) with the included 
matter field can be interpreted as the Jackiw- Teitelboim action (without a 
cosmological constant) in the non- commutative geometry of the "two sheet" 
manifolds Y x Z2, where Y is the 2-d space-time and Z2 is the internal space 
containing just two points |J§], |4"0|| . The matter field plays the geometrical 



role of the distance between the two points in the internal space. 

Our present model has also connections to the string theory on the curved 
backgrounds and to the exact 2-d conformal field theories. Indeed, in the 
conformal gauge, the classical action reads 

S = -— [ d 2 £(-d//<90 + K<pdQ j dQ 3 ). (113) 

This is obviously an action of the non-linear cx-model where /x, and Q- 7 are 
the coordinates of the target manifold with the metric 

ds 2 = -dfid(f) + K(j) dQ j dQ j . (114) 

It is not difficult to see that the metric ( |114|) (it is written in the so called 
Rosen coordinates) describes a single gravitational plane wave propagating on 
N + 2- dimensional target! |29|, |3(| In other words, a single gravitational wave 
in N + 2-dimensional target yields the cx-model action describing collisions 



of the two gravitational waves in two dimensions. Generalizing the work pi 
Brooks has shown that adding the target dilaton background in the critical 
target dimension to such a cx-model, one obtains an exact conformal field 
theory ||42j| . It would be interesting to study our present model from this 



point of view. All the mentioned features of the model look quite promising 
for further investigations and we shall certainly return to those problems 
elsewhere. 

I thank J. Frohlich, K. Gaw§dzki and R. Jackiw for enlightening com- 
ments. 



3 This comment is due to K.Gaw§dzki. 
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6 Appendix 



In this appendix, we evaluate the functional derivatives of the determinant 
in (pq), which were needed for obtaining the explicit form of the semiclassi- 
cal field equations. In evaluating the Traces we carefully keep in mind the 
definition of the scalar product ( f)0|) . 

Start with the derivative with respect to \i c . 



5 



-<e|exp{- E [-^e-^ia(^S)]}|? 



> . 



(115) 



We have used the heat-kernel regularization 44j based on the following 
representation 



POO (\ T ^ 

lnx = — / — e TX + (an x independent constant) + 0(ex) 

Je r 



(116) 



The "bras" and "kets" in Eq. (|115|) has to be understood in the standard 
sense. The asymptotic expression for the heat kernel for small e was obtained 



for an arbitrary elliptic operator in two dimensions [[44]] . If M has the form 

1 



M 



;V« + B a )^g^{sjp + Bp) - B { 



-9 



thenF] 



-9 + ^RV^+lB^ + Oie). 



117) 



;ii8) 



4 We put the sign + in front of R (see also Alvarez |4j| Eq. (4.38)), because R is given 
by Eq.©. 
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Performing the Wick rotation to the Minkowski time, we can use Eq. ( |118| ) 
for evaluating the heat kernel ( |1 1 5| ) . In our case 

B a = ^d a hi\<l> c \, (119) 
B = e~*>[ - \d 2 In |0 C | - i(01n |0 C |) 2 ], (120) 

hence 

^ Met i- 4 e " P i a( ^»l = i^ +I ^[2#.n|0 c | + (S ln|0 c |f]. 

(121) 

Note that the functional derivative with respect to fi c is the local expression. 
We also did not consider the first term in the r.h.s. of Eq. (118), which 



is eventually to be cancelled by the two-dimensional cosmological constant 
counter-term. 

Next we compute the variation with respect to <p c . First of all we note 
that Eq.([L2"lD implies 



lndet[-^e-^«9(0 c (9)] = 

1 I d 2 e(T^^V c +T^/i c [29Mn|0 c | + (ain|0 c |) 2 ]) + F(0 c ),(122) 

J 407T lD7T 

where F((p c ) is some /^-independent functional. But Eq. ( |122| ) itself gives 

l n det[-^d(0 c d)]=F(0 c ). (123) 
2h <p c 

Hence, 

5 , , r i „ 1 „, , i 



-lndet[ - — e -^—d(<j) c d)] = — t[&IIc ~ d(fi c d\n 



+ 5W) lndet ^ ~ ^ + (9ln l<Pcl)d] }' (124) 

Now we have 

51ndet{ - ^-[d 2 + (<91n|0 c |)<9]}| 0e=t = (125) 
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fOO Arp i 

-6 1 -TVexp { - r [ - + (d\n\M)d]]}U.., = 

f drTt { - i' 9 *) 3 ^ {^< a - '>< a < 2 + 7 s 4< a + »«} 

where a > is a (small) "Euclidean" cut-off, damping the oscillatory be- 
haviour of the exponent. Now we wish to evaluate the last Trace in Eq. (|125|) . 
Form the basis of the space of fields Q as follows 

#±,fc,pM = 6(±t)J (kt)-^=e- ip °, k > 0, (126) 

v 2ir 



where 6(t) is the usual step function. Using the theory of the Hankel trans- 
formations it is easy to establish the relations of orthogonality 

J r2 dtda t KkA^^k'At,*) = (-l)^^(A; - k')5(p-p') (127) 

/ dtdatilj* ±ktP (t,a)i> Tk/>p ,(t,a) = 0. (128) 
and the relation of completeness 

poo poo | 

/ kdk / d P (*; fcp (0 WO - *- kp (0*-k P (O) = -M - 1')5(* - a'). 

JO J—oo Z 

(129) 

(Note that in the r.h.s. of Eq. (|129|) stands the 5-function 5(£, £') with the 



respect to the inner product (|50"D for ji = and = t.) Therefore, the Trace 
of an operator O is given by 

p poo poo 

TrO = J R dtdat j o kdk J_Jp(** +kp (0O* +kp (0 - ^ kp (0O^- kp (0) 

(130) 

Using Eq. (|130|) , we can easily evaluate the Trace in Eq. (|125|) . We have 
51ndet{ - ^[d 2 + {d\n\<P c \)d]}U c=t = 

/ dr / dtdat kdk dpexp (- rf— (a - i)A: 2 + — (a + i)p 2 }} 
Je Jr 2 Jo J-oo 1 2n 2n > 
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roo r roo roo r T 1 

= / dr dtdat kdk dpexpj -[(a - i)k 2 + (a + i)p 2 ] \ 

Je JR 2 JO J-oo 1 2Ti > 

^(^)^(*w(o*wo - *vo*-*(0)- (131) 

Now the formula (see [p7| ) 

f°° 2 2 1 ( a2 +m 2 ) iAu 

/ rdre-" r J (Ar)J (pr) = —^^^ U^) (132) 
Jo 2p z 2p z 

and the standard Gaussian integration explicitly give the integrals over k and 
p. We obtain (for t ^ 0) 

51ndet{ - ^-[«9 2 + (9 In |0 c |)d]}k=t = jH dr dtdat^- 



-4, ft l ^ e-T^FJ ( J (fl 2 (t)-fl 2 (-t))l. (133) 

: + i)2r7T L(a-i)r v r(a-i) A w v v ; 

We can rewrite Eq.( |133| ) as follows 

^_lndet{-^^ + (.ln| 0c |).]}|^^ 

= -[9{t) - 9{-t)\ [' dpp- 1 ' 2 1 Uto(e-^Jo(^))' 

4y/2n(a - i) * L a - i J 

i 1 fit 2 

= ~W) - OH)] 1 \ -A^TT\ L S dpp- 1/ V^TJ (-^ T )) 

A^2ix(a - i) * KM a ~ 1 

= , 1 W r 2/£ dpp- 1 / 2 e-^J (^ T ). (134) 

4^/27r(a - i) * * Jo a-i 

Now we can decompose the integral over p 
r^ 2 /e dp _^ ip 

/ — e — « J ( :) = (135) 

jo y/p a — i 



rfrt 2 /e dp ip . a-i 

const + / — e a - i (J ( :) 7S= eQ " i 

Ji y/p a — l y/2np 



rht 2 /e (\p j a _ j 
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From the asymptotic behaviour of the Bessel functions |^7| , we conclude that 
the first integral in the r.h.s. of Eq.( [135| ) is convergent for e — > 0. Hence 

^ M e t {-^ + (a in W)a ]}|^ = 

i 1 1 fit 2 i ft 2 

- — -d t t dt- In (— —) = — In (136) 

8nt t V const ; 87rt 3 K eQ h K ' 

where Q is a finite constant. 
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